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unitarity and reversibility of microscopic dynamics and connect various statistical characteristics of 
many consecutive (continuous) observations in a quantum system subjected to external perturba- 
tions. Consequences of these relations in respect to the earlier suggested stochastic representation 
of interaction between two systems are considered. 

PACS numbers: 72.70.-|-m, 05.60.Gg 



I. INTRODUCTION 

Both the classical and quantum mechanics is time- 
reversible and unitary (conserves classical phase volume 
and quantum probability). In thermodynamic systems 
(ensembles), the unitarity manifests itself in strong con- 
nections between noise and dissipation, while the re- 
versibility in time symmetry of noise and reciprocity of 
transport processes. The famous examples are the Ein- 
stein relation [ij , the Nyquist formula , and the On- 
sager reciprocity relations . All these are the relations 
between (i) second-order (quadratic) correlators of equi- 
librium noise and (ii) linear parts of complete, possi- 
bly nonlinear, responses to external perturbations (the 
expansion of the responses in a series over powers of 
"perturbing forces" is meant). Later, the fluctuation- 
dissipation theorem (FDT) proved by Callen and Wel- 
ton Q and the Green-Kubo formulas 0, 0] exhausted 
this linear theory. 

The first nonlinear generalizations were obtained by 
Efremov who proved the quadratic FDT, which con- 
nects (i) quadratic parts of the responses and (ii) third- 
order (cubic) equilibrium correlators ((iii) linear re- 
sponses of quadratic correlators also take a part here, 
but can be excluded). The next, fourth-order, relations 
connect together (i) cubic components of the responses, 
(ii) equilibrium fourth-order correlators and, besides, (iii) 
quadratic responses of quadratic non-equilibrium corre- 
lators ((iv) linear responses of cubic correlators are in 
play too, but can be excluded). Their investigation was 
started by Stratonovich Q who found that "cubic FDT" 
does not exist (for details and more references see e.g. 
Ullllll)- Nevertheless, the fourth -order relations are 
much useful, for instance, when analysing low-frequency 
fluctuations in transport and relaxation rates, especially 
flicker fluctuations 

The producing formulas for the whole (infinite) chain 
of arbitrary-order fluctuation-dissipation and reciprocity 
relations (FDR) were obtained in j^. In the works 
|l2l ITsI l these results were extended to non-equilibrium 
steady states of open systems. In the exten- 
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sion to the thermic perturbations was developed, that is 
perturbations of a probability distribution (density ma- 
trix) of the system, in addition to perturbations of its 
Hamiltonian which are termed dynamic ones. Exam- 
ples of various a ppli cations of FDR can be found e.g. 

in El p Eiiii EE m III m 

It is desirable to combine all the infinite variety of 
arbitrary-order FDR into a compact visual generating 
FDR for the probabilistic functionals or corresponding 
characteristic functionals. In the framework of classical 
mechanics, this was realized in 

[HEl (for review, see 

In the quantum theory, time-differed values of any 
interesting variable X{t) (an operator in the Heisen- 
berg representation) do not commute one with another, 
X{ti)X{t2) ^ X{t2)X{ti). But its measurements in 
macroscopic devices are subjected to the classical descrip- 
tion language, being thought as a commutative stochas- 
tic process, x{t) , whose values are usual c-numbers. 
Hence, neither probabilistic nor characteristic functional 
of x{t) has a sense, until a concrete definition of all x(t)'s 
correlators (statistical moments), in terms of the X{t) , 
is chosen. 

In general, the two lowest-order correlators only seem 
unambiguously defined: 

{xit))=TTX{t)po, 
{x{h)x{t2)) = Tr {X{h) o X{t2)) po , ^ ' 

with Po being the statistical operator (density matrix), 
and o designating the symmetric product (Jordan prod- 
uct), AoB = {AB + BA)/2. The subscript "0" of po 
means that po represents a quantum statistical ensemble 
at a fixed time moment, e.g. t = 0. The right-hand sides 
present definitions of the angle brackets on the left sides, 
that is effective statistics of classical image, x{t) , of the 
quantum variable X(t). 

What is for the higher-order correlators 
{x{ti)x{t2)-.-x{tN)) , unfortunately, at any N > 2, 
there are N\/2 > 1 different symmetrized (Hcrmitian) 
expressions produced by various permutations of X{tk) 
in ^TTiX{h)Xit2)...XitN) + XitN)...Xit2)Xih))po, 
plus uncountably many their weighted linear combina- 
tions. 

Moreover, in fact even the second correlator can be 
introduced in an alternative way. For example, let us 
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define the characteristic functional (CF) of x{t) by the 
identity 



exp 
Tr exp 



P^v{t')X{t')dt' 



In Po 



(2) 



Here t > and v{t) is an "arbitrary probe function" 
(test function). Double differentiation of (0) by v{ti) and 
v{t2) at v{t) = gives 



x{ti)x{t2)) = Tr / X{ti)p^X{t2)pl-''da 



(3) 



and A'^-order moments include iV different pns ] 
whose sum equals to unit. As pointed out in [l^l, under 
this specific definition of the CF all the FDR between the 
angle brackets look absolutely similar in both classical 
and quantum case. 

Of course, to become practically useful, a choice of def- 
inition of the CF, i.e. definition of all statistical moments 
{x{ti) ... x{tj^)), should be based on analysis of real mea- 
surement procedures. The well known examples presents 
quantum theory of electromagnetic field fluctuations (see 
e.g. HSl)- During last decade, in the original works by 
Levitov, Lesovik, Nazarov and others [^[2^|2EI3| (see 
also references therein) the continuous measurements of 
the charge transport in electric devices were analyzed. 
The results of these works allow to conclude that fre- 
quently the adequate construction rule of the CF is the 
chronological symmetrized product: 

(exp [j^v{t')xit')dt'))^ 
Tr fecp (i v{t')X{t')dt'^ po ex^ (i /J v{t')X{t')dt'^ 

(4) 

Here fexp and exj5 are chronological and anti- 
chronological exponents, respectively. According to this 
rule (and to general properties of the trace operation 
Tr), the quadratic correlator remains as in For the 
higher-order correlators the rule Q prescribes 

{x{h)x{t2)x{t3)) = Tr iX{t,) o {X{t2) o Xm))) Po , 

{x{ti)x{t2)x(t3)x{t4)) = 

= Tr {X{h) o {X{t2) o (X(i3) o X{U)))) Po , 

(5) 

and so on, where, for definiteness, the inequalities ^at > 
... t2 >ti are presumed. 

Below, we will derive generating FDR for so built sta- 
tistical moments. All the more this is interesting be- 
cause the same rule (@J) naturally arose in the course of 
the so-called "stochastic re pre sentation of deterministic 
interactions" psl |23 | (the general method for 

correct construction of "Langevin equations" introduc- 
ing thermodynamic noise and dissipation into quantum 
or classical dynamics). 

The strong argument for the benefit of the rule Q is 
that it can be deduced from the correspondence principle. 
To see this, firstly consider classical systems. 



II. CLASSICAL CHARACTERISTIC 
FUNCTIONALS 

Let a classical system has canonical variables (coor- 
dinates and momenta) F = {q, p}. Generally, the sys- 
tem undergoes a dynamic perturbation from its outside, 
which means that its Hamiltonian, Ht = Ht{T) , is time 
dependent. Introduce also the corresponding Liouville 
operator Lt , and the evolution operator Zt : 



gxp 



Lrdr 



Eventually, we are interested in the evolution and fluctu- 
ations of variables X which represent definite functions 
of the phase space point currently occupied by the sys- 
tem: X — X{T) (breafly, "phase functions" ) . 

The essential properties of the Liouville and evolution 
operators are as follow: 



/ A{T)LtB{T) dT = -J B{T)LtA{T) dV 
J Ztp{T) dr = J p(r) dV , 

z^'x{T)Zt = x{z^'r) , 



(6) 



where ^(F), B{T), X{T) and p(r) are any phase func- 
tions (such that ^(r)i3(F) and p(F) are integrable), and 
Ft(F) stand for the current values of the canonic vari- 
ables (at time t ) expressed through their initial values 
F (at the initial time moment t = 0). In other words, 
Ft(F) is the solution of the Hamilton equations under 
initial condition ro(r) = F. 

Importantly, the latter equality in © is operator 
equahty, that is both A:(F) on the left and X{Z^^T) 
on the right have the sense of multiplication operators. 
Combination of this equality with the previous one im- 
plies 



Z^'X{T)Zt ^ Xit,T) = XiTtiV)) 



(7) 



Hence, the composite operator Z^^ X(T)Zt reduces to 
operator of multiplication by the time-dependent number 
X{t^Y) , which is nothing but trajectory of the variable 
X under initial conditions F. 

Next, let us be convinced that CF of any variable 
(phase function) X{t,T) = Ar(Ft(F)) can be expressed 
jlOj by the formula 



exp 



/o v{T)x{T)dT 

/^(/o*[L, + w(T)X(F)]dr) po{T)dT 



(8) 



Here po is the statistical operator, i.e. distribution func- 
tion, of the system at time t — Q, and again v{t) is arbi- 
trary probe function. In the angle brackets, x{t) means 
the internal variable X{t,T) as perceived by an outside 
observer (which knows nothing about F) and interpreted 
by him as a stochastic process. 
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To justify it is sufficient to make standard "dis- 
entangling" of the complex exponent in the lower row of 



/ fecp [j^[Lr + viT)X{T)] dr\ poiT)dT = 

J Zt^[j^v{T)z-^x{r)ZrdT] po{T)dr 

= / fecp [/o* v{t)X{t, r) dr] po{T)dT 



(9) 



Here, the latter transformation is made with taking into 
account the identity 121 and the second property from 
© ■ Evidently, final expression in Q coincides with what 
is meant in the angle brackets in ©, i.e. the CF of 
the path X{t,r) whose uncontrolled dependence on the 
initial conditions F turns it into a random process, x{t) . 

According to (jSJ-©, evaluation of the CF is equivalent 
to solution of definite differential equation: 



exp j'^ v{T)x{T)dT^ = J pdT , 
dp/dt - [Lt + vit)X{r)]p , 



(10) 



where the function p = p{t, F) satisfies the initial condi- 
tionp(0,F) =po(r). 

In principle, all what just was said is known. The 
representation ifTUI) . or (jSJ, used in liol lisll . is variation 
of so-called Feynman-Kac formulas which connect 
path integrals and differential equations. In fact, this 
representation of the CF is valid not for deterministic 
evolution only, but also for Marcovian stochastic evolu- 
tions. In this case, F designates instant state of a (multi- 
component) Marcovian random process, and Lt its evo- 
lution (kinetic) operator 
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III. QUANTUM CHARACTERISTIC 
FUNCTIONALS 

Following the correspondence principle, it seems natu- 
ral to suggest formulas © and/or l(Tn|l be the basis for 
definition of the CF in quantum case. 

Now, the Liouville operator changes to the commuta- 
tor: Lt^ = Ht]/h (with [A,B] = AB - BA, and 
^,A,B being arbitrary operators). Quantum analogue 
of the Zt exploited in previous section is super-operator 
whose action is defined by 

Ztp^U{t)pU-\t) , Uit) = loHrdr] 

The phase functions X{T) and X(Ft(F)) are replaced 
by the operators of quantum variable (observable) in the 
Shrodinger and Heisenberg representation, respectively, 
X and X{t), where X{t) = Z'^X = U-\t)XU{t) . 

What is for the operation of multiplication by X{r) 
in 1)10(1 . it should be replaced by super-operator of the 
symmetric product: X(F)$(F) ^ X o ^ . The matter 
is that it is very hard to suggest something else. Under 
this extension, the equality ((TJ remains valid: 

z^r^ix o (Zt$)) = x{t) o $ = {u-^{t)xu{t)) o $ 



(but now, generally, the super-operator X{t)o in none 
base reduces to scalar multiplication). 

Thus one comes to the construction of quantum CF as 
follows: 



yexp 



dp 
'dt 



v(T)x(T)dT 



Trp , 



^-[p,Ht]+v{t)Xop 



(11) 



(12) 



with the initial condition p{t — 0) ~ po- One can easy 
verify that substitution of the formal direct solution of 
into produces just the formula Q, with X{t) 
standing for Heisenbergian operator of the quantum vari- 
able and x{t) its effective commutative image. 

The doubtless advantage of such definition. ([TT |) - ((T^ . 
of quantum CF is its differential nature |25l l2Qj , which 
highlights its automatic agreement with the causality 
principle too. The formulas © demonstrate that be- 
cause of the causality all the corresponding higher-order 
statistical moments appear asymmetric with respect to 
time inversion. It is natural: earlier observations (mea- 
surements) can affect later ones, but opposite influence 
is impossible. 

Further, we want to extend to an arbitrary set 
of variables under simultaneous continuous observation. 
With this purpose, let us replace v{t)X in H12|) with the 
operator 



Vt 



'X„ 



x„ 



(13) 



where \p) are states, or vectors (in the Dirac's desig- 
nations) , which constitute a complete orthonormal base. 



and 



1 1'^ 



arbitrary probe functions. This is most gen- 



eral form of the observation. Analogously, most general 
external perturbation can be described as 



Ht — Ho ~ Ft , -Pt ~ ^ ft^^tiv , 



(14) 



where Hq is Hamiltonian of the "free" system, and ft 
are "perturbing forces" . 

Then, instead of p2|l and we have 



(15) 



exp 



Tr p 



S(F,; Ho.Fr) 
(16) 

From here, the repeated indices p or v mean summation 
over them; x^i,{t) are effective classical images of quan- 
tum variables X^^{t) = U-\t)X^^U{t) ; E{Vr;Ho, Fr) 
will be used as shortened designation of the CF, i.e. the 
angle brackets. 

The second and third arguments of S, as well as 
the subscript under angle brackets, remind about eigen 
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Hamiltonian of the system and its perturbations. Direct 
solution of lfT3|) yields 



= Tt t^i^ /J [-j-{Ho - Fr) + \Vr] drj poX 

XCX$lj^[j-{Ho-Fr) + ^Vr]dT) = 



(17) 



Tr fexp 



exg i/owrX^,(T)dT 



The S's arguments Vt and Ft can be understood as 
either two sets, of the probe functions v^'^ and the forces 
Z/^" , or the whole operators introduced in (|13|l and (|14|l . 
Of course, the pair x^^it) and x^^{t) in H16|l can be 
interpreted as two mutually conjugated complex-valued 
random processes (at least while Ft assumed Hermitian) . 



IV. THE GENERATING FDR 

In the present paper we confine ourselves by the as- 
sumption that the initial density matrix, po , is canonic 
thermodynamically equilibrium one: po oc exjp^—Ho/T) 
(in the classical theory, not only this case already was 
considered [TH ITsj l but also the case of non-equilibrium 
initial distributions 0, IT^). 

Following the recipes of [ij, let us make three trans- 
formations of the operator expression placed after the 
trace symbol in H17|) . 

(i) Firstly, transpose this expression as a whole, with 
the help of the usual rule {AB...CY = C'...B'A' , where 
the prime means the transposition: A' ~ (A^)* = (A*)^ 
(symbols f and * will stand for the Hermitian and com- 
plex conjugation, respectively). As is well known, this 
operation does not change the trace. 

(ii) Secondly, invert the direction of the time account, 
by means of rewriting chronological exponents as anti- 
chronological ones, and vice versa, as in the examples 



exr 



IoA{T)d, 



6xp 



IoMt^r)d7 



(iii) Thirdly, rewrite po in the form ^/po^/po and 
"drag" one of these multiplicands to the left and another 
to the right and after that again unify them, with the 
help of the trace property Ti ABC = Tr CAB . 

The result of these manipulations is an expression 
which looks quite identical to the initial one, but with 
modified operator of the observation instead of Vt and 
modified operator of the perturbation instead of Ft . To 
write the result, of course, it is convenient to choose the 
states 1^) as a complete set of orthonormal eigenvectors 
of the free Hamiltonian Hq : Hq\p) — E^\ij), with the 
eigenvalues (energies) E^. Besides, introduce two oper- 
ators 

Lo$ = Ho\lh , Uo{t) EE exp[-zi/ot/7j] , 



the matrices (in terms of the chosen basis) 

2T 



S„f, , C^f, = sinh , cosh , A^^ = ^ tanh ^ 



(18) 

and super-operators C, S and A whose action is de- 
termined by these matrices: 



C = cos i^Lo) , (C$)^^ = Cf,^<^^^ , 
'^Lo) , {S<f),,,= 
(A<I>)^^ = A^^$^^ 



S - sin {^Lo) , (S$),,, = 5,,,$^, , (19) 



Then, finally, we can formulate the generating FDR. 
It is expressed by the formulas as follow: 

E{Vlr ; K FU) = 2(K ; Ho, Fr) , (20) 

whereby; = vt^'^X'^, , Vr = v'^'^X^, , F^ = f^X'^, , 
Fr = fl^'^X^,, , and 



2-q r 



or, in another equivalent form. 



" fr ' 









' Fr' 




■ c -^s' 




' Fr' 


.Vr. 




.is c . 







(21) 



(22) 



For convenience, the observation and perturbation vari- 
ables are unified into the column vector. 

The combined transformation ^r ^ ^'t-r represents 
the time reversal. Thus, (^01) and if^ or (1^ state in- 
variance of the CF under (i) simultaneous time reversal 
of both the probe functions and external forces and (ii) 
their mutual mixing as described by H21|) and (|22|l . 

Importantly, when dealing with H2U|) - (|21|I one should 
remember about the time translational invariance: any 
joint temporal shift of Vt and Ft does not change 
'E.{Vr] Ho, FrYa value, since po is invariant with re- 
spect to free (unperturbed and unobserved) evolution 
([po,i^o] = 0). 

It should be underlined also that by the very definition 
(see p3f) ') of the operators X^^ we have 



H'o\p*) ^ E^\p*) 



(23) 



where are eigenvectors of the transposed free Hamil- 
tonian Hq , with the same eigenvalues E"^ as that of \p) 
(one can write also \p*) — \p)* = {lA')- Hence, during 
the time reversed evolution the operator X'^^, represents 
those quantum variable which is represented by opera- 
tor X^^ in the direct process. Consequently, in terms of 
their effective classical images, the left side of (I2U|I looks 
as 

S{Vl_r;H'Q,F^_^) 



(e.vUlv7.r^,.{r)dT\) (24) 
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In other words, under the time reversal the matrix {v!^'^} 
behaves hke the operator Vr , that is it undergoes trans- 
position: ^ = Vt-T (similarly, f!^" ^ 

In the classical limit (formally ?i ^ 0), we have C ^ 1, 
S ^0, 2S/h -> T-^Lo , and formula ^ takes the form 



Fr^Fr , Vr= K 



T-^L^Fr 



(25) 



At that, the operators Fr and Vr turn into phase func- 
tions, Lq {y qH^ys/ p—{V pHq)'S/ q , and the prime in 
(|20|l means replacement of {q,p\ with {q,—p\. 

Comparison between (|22|l and H25() reveals that the ob- 
servations anyway are influenced by the perturbations, 
but opposite effects exist in the quantum theory only. 
For illustration, let us choose = , that is the real 
direct-time observations (described by Vr) are made in 
equilibrium system. Then, according to (pn| - lf^ . detec- 
tion of exactly time-reversed results (described by V/_^ ), 
in equilibrium system too, is in certain sense the same as 
detection of the direct results but under specific nonzero 
perturbations. This unpleasant peculiarity of the quan- 
tum theory is not surprising: since any measurement in- 
fluences subsequent ones, but not vice versa, the mere 
rearrangement of their results, generally speaking, could 
not realize under same conditions. 

Various particular FDR can be obtained from l|2U|l - H22|) 
by either some special choice of Vr and Fr or varia- 
tional differentiations with respect to v!^'^ and f!^'^. Ev- 
idently, if both Fr and Vr are Hermitian, then both 
transforms Vr and Fr also are Hermitian. But if Vr 
is quite arbitrary, then one should be ready to deal with 
non-Hermitian perturbation Fr- 

For example, let us choose Vt = 0, i.e. there is no 
observation at all. According to l(T^ or ((TTIl . of course, 
in absent of observations the evolution of the statisti- 
cal operator p is unitary under whatever perturbations. 
Therefore Trp = Trpo = 1 , and S(0;iJo, ^V) = 1- 
Consequently, Ip^) and if^ . together with ifT^ . yield 
the equality 



1 = S (^^ALoFr ; Ho, Fr 
or, in other equivalent designations. 



(26) 



exp 



J|tanh(^|^)/rx,.(r)dr 



(27) 

These relations are valid at arbitrary perturbation oper- 
ator Fr . The subscript under the angle bracket reminds 
about the perturbation. 

To make Eas l26l27l better transparent, combine them 
with the identities 



LoFt — j-[Ft,HQ] 



-[Ho-Ft,Ho] = -Pt 



where Pt is operator of the instant energy power being 
pumped by the perturbations (notice that dHo{t)/dt = 



d{U-^{t)HQU{t))/dt : 
rewrite Eq|2S|also as 



U-'^{t)PtU{t)). Hence, we can 



-ffo, Fr] 



(28) 



with arbitrary Fr and related Pt — j^\HQ,Ft]. 

In the classical limit, A disappears, in the sense 
that A^i, 1. Therefore, (^5)1 turns into equality 
(exp(-A/T)) = 1 where A is total work pro- 

duced by perturbations during time interval (0, t) . More 
generally, perhaps, APt/T can be interpreted as opera- 
tor of entropy production. 



V. PROBABILITY FUNCTIONALS 

Provided statistical moments of quantum variables 
X{t) and their CF are defined (by H15|) - (|17|l . in our case), 
we can introduce the probability functional (PF) of their 
classical (commutative) equivalents x{t) . We will des- 
ignate it by W{x{t) ; Hq, Fr) ■ As usually, is represents 
functional Fourier transform of the CF: 



W{x{t) ■ Ho, Fr) = / exp -i <-a;^,(T)d- 



xEitVr;Ho, Fr)dVr 



(29) 



where dVr = Ylr / '^^) ■ 

Let us apply this transform to both sides of the 
FDR (pn|) . again with the indices being related to the 
eigenstates of Hq. We omit rather bulky manipulations 
and write their result in two steps: 



W{x\t-T);H'„F^_,) 



(30) 



exp 







— / — tanh 


. Jo ^ ^ 





X W{xiT) ■,Ho,Fr) 



Wix{T); Ho, Fr) 



(31) 



exp 



h 
2i 



dr^^S^.C^, 5/r 5x 



(r) 



X W{C-'x{t) ■ Ho, C-'Fr) , 

where (C^^*)^^ = ^^,u/C^,^. 

In these two formulas, both the exponents in mid- 
dle rows result from the left bottom and right top non- 
diagonal elements of matrix (or (j^ ). respectively, 
i.e. from mutual mixing of perturbation and observation. 
Thus, the exponent in (|31|l reflects disturbing action of 
observations, which is quite unpleasant peculiarity of the 
quantum case. Of course, in fact this operator-valued 
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exponent acts as an integral operator. We leave it un- 
wrapped till possible separate work (though example of 
similar operators can be found in [2a|)- For the present, 
confine ourselves by the sad conclusion that generally 
W{x{t); Hq, Ft) , and thus PF of the reversed process, 
W{x'{t—T); Hq, F[_^) , relates to PF of the direct process 
W{x{t)] Hq, Ft) in some non-local way, with respect to 
both /(^'' and Xfii,{t). 

However, it is not hard to notice that a degree of the 
non-locality is proportional to tanh^(£'^^/2r) , hence, 
it is negligible in respect to low-energy quantum transi- 
tions. In the classical limit {h 0) the exponent in l|31() 
turns into unit, W W , and reduce to the 

purely local relation 



W{x'{t- 



eM-A/T)W{x{T);Ho,Fr) , 

(32) 



where A is again the work of the external forces /j 
and ^'{q,p) = i'fq, —p). This FDR is equivalent to what 
was obtained in Curiously, if quantum CF is de- 

fine d by E qE) instead of Eas ll5ll7l then EqES replaces 
Eas l^UlJTl even in general quantum case. 



VI. STOCHASTIC REPRESENTATION OF 
RESPONSE TO PERTURBATIONS 

The above consideration demonstrated a lot of formal 
symmetry between perturbations and observations. This 
symmetry suggests that the perturbing forces f^"^ can 
be treated as another test (probe) functions which cor- 
respond to watching for additional ghost variables. Let 
the latter be named as ?/^^ {t) . To define them, following 
El El 113, we rewrite Eas ll5ll7l in the form 



exp [vtf'-x^.iT) + fry, At)] dr) 

Tr fecp ( /J [-j-{Ho - Fr) + \Vr\ dr) poX 
x^(/o [i(^o-^^.) + iK] dr) 



(33) 



Essentially, it is assumed here that an imaginative prob- 
ability measure hidden behind the angle brackets is itself 
independent of the forces. In other words, all the random 
processes x^^it) and y^v{t) are meant be characteristics 
of the unperturbed dynamics governed by the Hamilto- 
nian Hq. The subscript "o" serves to remind of this 
circumstance. 

Therefore, it may be convenient to rewrite H33|l once 
more, in the form 

(exp /„* [v^x^^ir) + frVf^Ar)] dr) = 
= E{Vr;H„,Fr) = 
= Tr ^ (/J [j-fr + i<1 X",uir)dT) poX ^^^^ 

X ^ (/o [-l/r + 5<1 X",uir)dT' 



where X°^{t) are the operators X^^ considered in the 
interaction representation, thus representing free evolu- 
tion: 

According to (O-JMI, if V"f ee then 

Jexp /g* fry^.Ar) dr) ^ = S(0; H^, F,) = 1 , (35) 

that is any statistical moment of y's themselves is equal 
to zero. But their correlations with x's differ from zero: 



J|a;(tj)?/(7 



n 



<5/(t„ 



Fr=0 

(36) 

(indices p.,!/ are omitted). Hence, cross-correlation be- 
tween N copies of X and M copies of y represents M- 
order response to perturbations of an A^-order statistical 
moment of the x's. 

Interestingly, the relations (|35|l are valid also for the 
x's and ?/'s cumulants (semiinvariants) whose generating 
function is In E{Vt; Hq, Fr). The proof is trivial: 

In (exp / [vrx{T) + frVir)] dr) ^ = 
= [eicp{JdTU5/5gA In S(K; i?o,G,)]G^^o 

(again without indices) . Consequently, instead of H3t)|) we 
can write 



]^a;(ij)?/(7 



n 



5f{T,n) 



\{x{t, 



(37) 

with the superscript "c" marking the cumulants. 

The union of the two sets of random processes, quite re- 
alistic x's and indeed rather illusive y's, gives stochastic 
representation of the system's response to perturbations. 
If the latter are caused by interactions with some other 
dynamical system "D" , then we make first step towards 
the stochastic representation of deterministic (quantum 
or classical) interactions, which was suggested in E3 and 
developed in E1E3,E1|- For instance, our system can 
serve as thermostat for "D" . General stochastic equations 
which describe "D" under influence by the thermostat in- 
evitably include the y's whose main effect is dissipation. 



VII. TIME REVERSAL AND GENERATING 
FDR IN THE STOCHASTIC REPRESENTATION 



Combining Eas l20l2^ with Eg 1331 one can simply re- 
formulate the generating FDR in terms of a;'s and ?/'s : 



(38) 



exp /J [v^x'^At -r)+ fry'^^At - r)] dr) ^ = 
= (exp ll br5?^.(r) + !ry,Ar)\ dr 



7 



where, of course, x'^^ = Xi,^ , = y^^i, and x, y relate 
to X, y absolutely similar to (|2 l|l : 



2 '-'a"^ 



(39) 



Evidently, it the matrices x — {Xfj_^} and y — {y^^i/} are 



Hermitian, x 



t - 



and — y , then the transformation 



(|39|l does not damage this their property. 

Taking into account po's invariance with respect to 
free evolution, the FDR H38|l can be expressed also in the 
form 

Xi^fj^ih ~t):< Xf,^{t) , y,y^{tQ - i) X y^^(t) , (40) 

where to is arbitrary time shift, and symbol x means 
statistical equivalence of left- and right-handed random 
processes (i.e. equivalence in the sense of statistical mo- 
ments, {x'{to — ti)...x'{t() — In)) ~ {x{ti)...x{t]y)) , and 
so on). 

Following [2^, it may be convenient to introduce 
another random processes, whose matrices are non- 
Hermitian: 



at) EE X{t) + f 2/(0 , 



(41) 



ri{t)^e{t)=x{t)-%y{t) 
In their terms the generating FDR look most simple: 

C.M(io-t) X exp(S^,/2T)ep.(t) , (42) 



(43) 



To some extent, the ^'s and 77's can be thought like 
amplitudes of quantum jumps, hence, their squares, 
and |?7p, like corresponding probabilities. Then, formu- 
las (|42ll and (|43|) reduce to familiar relations between 
probabilities of mutually reversed jumps. 

It is necessary to remember, of course, that generally 
(when Hq Hq) the left and right-hand processes in 
(|4()|l . (|42|l and 14,3|l relate to somehow different systems. 

Combining these statistical equalities with H36() (or 
(127))) and (or (|34|l '). and besides with the causal- 
ity principle, as well as with independence of statistical 
moments on to , one can construct relatively simple algo- 
rithms for derivation of many particular FDR. However, 
that are tasks for separate work. 



VIII. CONCLUSION 

To resume, we obtained generating fluctuation- 
dissipation relations (FDR) for a dynamically perturbed 
quantum system, assuming special but theoretically and 
practically important definition of its quantum statisti- 
cal moments (or corresponding characteristic functional) 
which describe consecutive or continuous measurements 
of the system. 

In addition, short and expressive formulation of the 
FDR in terms of the stochastic representation of quantum 
interactions was done. 
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